Stability of quasi-two dimensional zigzag carbon and its reaction pathway to graphene 
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Using the density functional theory we investigate a quasi-two dimensional carbon allotrope, ZzC, 
formed by square carbon lattice buckled in zigzag manner. By analyzing the Kohn-Sham energy and 
phonon dispersion obtained by lattice dynamical calculations we show that ZzC is stable with binding 
energy of 7.46 eV per atom. To examine the possible route of formation we find out reaction pathway 
from ZzC to graphene using nudge elastic band method, generalized for solid state calculations. The 
reaction pathway shows the formation of carbyne as the intermediate state. Such a pathway is seen 
to exhibit two transitions states with reaction barriers of 0.21 eV/atom from ZzC to carbyne, and of 
1.19 eV/atom from graphene to carbyne. Although ZzC is stable, upon hydrogenation it dissociates 
and prefers the carbyne-like structure by forming chains of polyacetylene. 

PACS numbers: 61.48.Gh,81.05.Fb,81.05.ue,88.30.R-,73.22.Pr,71.15.Mb,71.20.-b 



I. INTRODUCTION 

Variety of sp hybridizations of carbon allows it to have 
diversity of structures in all three dimensions 1 . After 
discovery of graphene, 2D family of carbon allotropes 
became subject of special interest. For example, hydro- 
genated graphene, called as graphane is under investi- 
gation by both theoretical and experimental groups^—. 
Other materials like graphone and graphYn have also 
been predicted theoretically. Zhou et at 10 showed that 
semihydrogenated graphone can act as a ferromagnetic 
material. On the other hand graphyne has been claimed 
to be better than grapheneii — . An exhaustive theoret- 
ical study was carried out recently by Wen et af who in- 
vestigated various allotropes of Group 14 elements using 
Density functional theory (DFT). They studied a variety 
of ID, 2D and 3D structures of group IV elements; par- 
ticularly their geometries and stabilities by means of for- 
mation energies. Their extensive study shed the light on 
various structures and their relative stability. They also 
found that in 2D a wavy square sheet of carbon is consid- 
erably more stable than other crystals except graphene. 
It thus becomes very intriguing to investigate such wavy 
material in details, which we call Zigzag carbon (ZzC). 

In this work we investigate the crystal structure, stabil- 
ity and electronic structure of ZzC. The relative stability 
is investigated by using binding energy (BE) and phonon 
dispersion spectra as BE cannot be considered conclu- 
sive. Moreover, we aslo investigate the reaction path- 
way from ZzC to graphene. Such a transition involves 
both atomic and cell degrees of freedom. We address 
this transition using Linear Synchronized Transit (LST), 
Quadratic Synchronized Transit (QST)i^ and General- 
ized Solid-State Nudge Elastic Band (G-SSNEB) 15 meth- 
ods. Although LST is easy to implement, QST, and spe- 
cially G-SSNEB are more accurate in predicting transi- 
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tion state. We observed that QST reaction barrier from 
ZzC to carbyne is in 5% in agreement with G-SSNEB 
one. 



II. COMPUTATIONAL DETAILS 

We perform constraint free full unit cell optimization 
using Generalized Gradient Approximation (GGA) with 
exchange correlation functional of Perdew, Burke and 
Ernzerho£i£ as implemented in Quantum ESPRESSO 17 
The primitive cell contains two carbon atoms. The unit 
cell optimization was terminated upon reaching the pres- 
sure cut-off of 0.1 Kbar. Kinetic energy cutoff on plane 
waves was set to as high as 1360 eV. The energy criterion 
for electronic self consistency was set to 10 -8 eV while 
that for structural optimization was 10~ 7 eV. For the 
optimization the /c-mesh was sampled using 16 x 16 x 1 
points using Monkhorst-Pack scheme, while the band 
structure was plotted on lines joining four symmetry 
points of Wigner-Seitz cell, and the individual line seg- 
ments were sampled using fifty /c-points. A vacuum space 
of 25 A was kept normal to the plane to avoid any inter- 
planar interactions. 

Phonon dispersion is obtained by lattice dynamical cal- 
culations performed using self-consistent density func- 
tional perturbation theory within linear response ap- 
proach. The energy threshold value for convergence was 
1.4 x 10 -15 eV. Force constant matrices are obtained on 
a 3 x 3 x 1 /c-point mesh of irreducible part of Brillouin 
zone, and are interpolated at arbitrary wave vectors. The 
dispersion spectrum is then plotted along the lines join- 
ing symmetry points sampled using total of 150 q points. 

For finding solid-state transition pathway between 
graphene and ZzC we used three methods: LST, QST 
and G-SSNEB. In LST all degrees of freedom (DOF), 
namely cell parameters and atomic positions, are var- 
ied linearly from react ants to products. The maximum 
on the resulting pathways corresponds to a transition 
state (TS). Better results of TS are obtained by QST, 
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in which one should first optimize the LST maximum 
orthogonally to its pathway in combined phase space of 
coordinates and cell vectors. (In this work we finally 
changed orthogonal optimization to unrestricted opti- 
mization because the resulting structure appeared to be 
another carbon allotrope with energy below ZzC.) Then 
a quadratic transit pathway is built from reactants to 
products through this TS. For detailed discussion on 
QST we refer reader to reference 14. The G-SSNEB 
method, which we implemented in QE, was developed 
by Henkelmann's groupie. The only difference with orig- 
inal method is the cell gradient (Eq.J2|)), which we used 
in the form given by Caspersen and Carter Choice of 
gradient is based on ease of numerical implementation 
and although the nature of gradients are different in both 
the cases; we believe that they should produce identical 
results in a simple case that we are studying. 

The gradient in the combined phase space is the gen- 
eralized forces, which are 



"^ = -"((- + W 1 ) T (2) 

where x, F, h, a are matrices of atomic positions 
in crystal coordinates, forces, cell vectors and internal 
Cauchy stress tensor accordingly. P is the external pres- 
sure, equaled zero in our calculations. 

III. RESULTS AND DISCUSSION 

We shall begin the discussion by the presenting the 
lattice structure and the stability, and proceed to elec- 
tronic structure and the reaction pathway from ZzC to 
graphene. We shall also talk about the effect of hydro- 
genation. 

A. Crystal: structure and stability 

Figure[T](a) shows the optimized crystal structure along 
with its lattice parameters and bond lengths. For sake of 
clarity various views are depicted, which also make the 
name zigzag carbon evident. The unit cell belongs to Cih 
point group of monoclinic crystal, with lattice parame- 
ters of a = 2.46 A, b = 1.65 A. As seen from the top 
view the bond length along the zigzag direction (1.51 A) 
is smaller than that along non-zigzag direction (1.65 A). 
These numbers are consistent with Wen et a? (1.50 and 
1.63 A respectively; bearing in mind the choice of dif- 
ferent treatment on exchange correlation functional). At 
this point we recall that the C-C bond length is 1.42 A 
and 1.62 A in graphene and SqC (planar square carbon 
allotropei^) respectively. It is thus intuitive to think that 
binding energy of ZzC may lie in between the two. In- 
deed, ZzC has the binding energy of 7.46 eV per atom 



which is higher than SqC (6.7 eV per atom) but not as 
high as graphene (9.6 eV per atom). However binding en- 
ergy alone cannot serve as the stability criterion. The sta- 
bility of ZzC becomes evident conclusively by analysing 
the phonon dispersion spectrum (Figure [U(b)). Obvi- 
ously, the lack of any imaginary phonons in entire Bril- 
louin zone strongly underlines the stability of ZzC. The 
phonon spectrum has total of six branches, three acoustic 
and three optical (out-of-plane, in-plane transverse and 
in-plane longitudinal for both types). It is interesting to 
note that although there is no band gap between acoustic 
and optical branches, an indirect band gap of 18.9 cm -1 
exists between the out-of-plane optical mode and in-plane 
transverse optical mode. Thus, from phonon dispersion 
it becomes evident that ZzC is indeed a stable material 
with a high binding energy. 



B. Electronic structure 

After validating the stability we now discuss the elec- 
tronic structure of ZzC. As seen from the band structure 
(Figure EJ^a)) ZzC shows a small band gap of 0.39 eV 
at the X' point of Brillouin zone. The band gap at V 
point is rather high (7.8 eV). We also note that ZzC 
does not show any magnetic character. Figure [2](b) also 
shows the charge densities of highest occupied and low- 
est unoccupied bands. Highest occupied band which is 
primarily p-type shows the stacking of baguette-shaped 
charge density, while unoccupied band completely embed 
the atoms within individual planes. It is interesting to 
note that both states show the derealization perpendic- 
ular to the zigzag direction. Looking at the delocalized 
nature of bands and bearing in mind that DFT tends to 
underestimate the band gap, one may speculate that the 
ZzC may act as a semiconductor with peculiar conduc- 
tance perpendicular to zigzag direction. 



C. Reaction pathway 

LST, QST and G-SSNEB pathways along with initial, 
transition, intermediate and final states, are presented 
in Figure 3. LST method does not have the interme- 
diate state found by QST and G-SSNEB independetly. 
It is explained by non-iterative nature of LST: it rep- 
resents only the shortest pathway between reactant and 
product, and is useful only as a quick and rough ap- 
proximation. Presence of the intermediate state, carbyne 
(weakly binded carbon chains), means a multistep reac- 
tion: along the pathway to the more stable graphene 
allotrope ZzC dissociates to carbyne. Carbyne structure 
(Figure 3) is composed of linear carbon chains^ weakly 
binded by dispersion interaction. Our results are in good 
agreement with parameters of betta-carbyne (or polycu- 
mulune type, cumulated double bonds) in Heinmann's 
work^i. In particular double C=C bond is equal 1.28 A 
in both cases, which is in agreement with the Steichen's 
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FIG. 1. (a) Optimized structure of ZzC. The name zigzag is evident from the side view. The dotted blue line indicate the 
primitive unit cell which has the parameters a = 2.46 A, b = 1.65 A and 7 = 90°. The primitive cell has two carbon atoms 
with the bond length of 1.51 A. (b) The phonon dispersion relation is shown along the symmetry points of Brillouin zone. The 
absence of soft modes confirms the stability. 
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(a) (b) 

FIG. 2. (a)Band structure and the DOS of ZzC. DOS are resolved as per their angular momentum components (s and p) and 
are plotted using the electronic smearing of 0.1 eV for better visualization. A small band gap of 0.39 eV is seen at the X 1 
point, (b) Charge density of highest occupied and lowest unoccupied bands show the declocalization perpendicular to zigzag 
direction. 



pricinple 22 . Interchain distances are 3.7 and 2.9 A m 
our and Heinmann's work accordingly. The bigger value 
in our case is due to different unit cells in the works. 
In our caclulation we used artificially small unit cell for 
carbyne: there are only two carbon atoms, whereas in 
Heinmann's work n-atoms cells (n » 2) are presented. 
Bigger cell allows more energetically favourable packing, 
which leads to higher interaction, thus smaller interchain 
distance. That difference also explains different transi- 
tion temperatures from graphene to carbyne: in our cal- 
culation reaction barrier from graphene to carbyne is 1.19 
eV/atom, that corresponds to about 4600 K (the method 
of temperature calculation is described in the next para- 
graph), whereas in the work of Nelson etc. 23 graphene 
was heated above 2600 K. The planar structure of car- 



byne is most likely an artefact caused by construction of 
primitive cell (to consider possibility of non-planar struc- 
ture supercell in the orthogonal direction to the chains 
should be constructed). In reality final structure could 
be a layer or even clusters of amorphous carbon com- 
posed of crosslinked linear chai mj 2Q i 24 . But this part of 
modeling lies beyond the scope of this article. 

To answer the question of stability of ZzC structure 
we calculated transition temperature corresponding to 
reaction barrier using rough approximation of atoms as 
independent simple harmonic oscillators having energy 
3kT each. Thus, a rough value of the transition temper- 
ature is E/(3k), where E and k are energy of reaction 
barrier per atom and Boltzmann constant accordingly. 
For transition from ZzC to carbyne it corresponds to 
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FIG. 3. Reaction pathways from ZzC to graphene using LST, QST and G-SSNEB methods (in the last one we used two 
climbing images because of two reaction barriers). 



around 815 K. Thus ZzC structure is a pretty stable 
composition. 



D. Hydrogenation 

At this stage we wish to point out an interesting ob- 
servation. Unlike graphene or SqG^, ZzC responds to 
hydrogenation in qualitatively different manner. We re- 
call that upon hydrogeation graphene forms grapheme, 
while SqC accommodates hydrogens by forming hexaco- 
ordinated bonds. In contrast, we observed that upon hy- 
drogenation, ZzC breaks down and forms stacks of poly- 
acetylene (-[CH] n -). The individual chains of polymers 
were spatially separated by more than 3.5 A, ruling out 
the van der Waals interactions. We verified the absence 
of any such interactions by enabling a correction term to 
the exchange correlation functional, to take into account 
the long range the dispersion i 25 i 26 Yet, such different re- 
sponse to hydrogenation is understandable by analyzing 
the nature of electronic structure. While graphene has 
one double bond, SqC has electron deficient bonding. On 



the other hand, all bonds in ZzC are single bonds and are 
completely saturated. By disintegrating into polyacety- 
lene carbon atoms gain significant binding energy (from 
7.46 eV to 8.31 eV per atom). Thus although the ZzC 
is extremely stable, we speculate that it can be easily 
converted into stacks of polyacetylene by hydrogenation. 
IV. CONCLUSIONS 

We performed the full unit cell optimization on zigzag 
carbon and found out that the ZzC may act as the semi- 
conductor with valance and conductance bands pecu- 
liarly delocalized only perpendicular to the zigzag direc- 
tion. Our calculations revealed that with binding en- 
ergy of 7.46 eV per atom ZzC is stable, which we also 
verified using lattice dynamical calculations. We also 
observed that ZzC breaks into chains of polyacetylene 
upon hydrogenation. Finally our reaction path from 
ZzC to graphene using QST and G-SSNEB showed that 
ZzC transforms into carbyne with reaction barrier 0.21 
eV/atom, roughly corresponding 815 K transtion tem- 
perature. We believe many other exotic 2D carbon al- 
lotropes have the same fate: breaking apart to carbon 
chains with following amorphous phase formation. 



1 Xiao-Dong Wen, ThomasJ. Cahill, and Roald Hoffmann. 
Chemistry A European Journal, 16(22) :6555-6566, 2010. 

2 D. C. Elias, R. R. Nair, T. M. G. Mohiuddin, S. V. 
Morozov, P. Blake, M. P. Halsall, A. C. Ferrari, D. W. 
Boukhvalov, M. I. Katsnelson, A. K. Geim, and K. S. 
Novoselov. Science, 323(5914) :610-613, 2009. 



3 R. Balog, B. J0rgensen, L. Nilsson, M. Andersen, 
E. Rienks, M. Bianchi, M. Fanetti, E. Laegsgaard, 

A. Baraldi, S. Lizzit, Z. Sljivancanin, F. Besen- 
bacher, B. Hammer, T. G. Pedersen, P. Hofmann, and 
L. Hornekaer. Nature Mater., 9(4):315-319, 2010. 

4 J. O. Sofo, A. S. Chaudhari, and G. D. Barber. Phys. Rev. 

B. , 75:153401, 2007. 



5 



5 D. W. Boukhvalov, M. I. Katsnelson, and A. I. Liechten- 
stein. Phys. Rev. B, 77(3):035427, 2008. 

6 Prachi Chandrachud, Bhalchandra S Pujari, Soumyajyoti 
Haldar, Biplab Sanyal, and D G Kanhere. J. Phys.: Con- 
dens. Matter, 22(46) :465502, 2010. 

7 B. S. Pujari and D. G. Kanhere. J. Phys. Chem. C, 
113(50):21063-21067, 2009. 

8 A. S. Barnard and I. K. Snook. J. Mater. Chem., 20:10459- 
10464, 2010. 

9 H. §ahin, C. Ataca, and S. Ciraci. Phys. Rev. B, 
81(20):205417, 2010. 

10 J Zhou, Q Wang, Q Sun, X S Chen, Y Kawazoe, and 
P Jena. Nano letters, 9(ll):3867-70, 2009. 

11 R. H. Baughman, H. Eckhardt, and M. Kertesz. J. Chem. 
Phys., 87:6687, 1987. 

12 A. N. Enyashin and A. L. Ivanovskii. Phys. Stat. Sol. B, 
248:1879, 2011. 

13 Daniel Malko, Christian Neiss, Francesc Vines, and An- 
dreas Gorling. Phys. Rev. Lett, 108:086804, 2012. 

14 T.A. Halgren and W.N. Lipscomb. The synchronous- 
transit method for determining reaction pathways and lo- 
cating molecular transition states. Chemical Physics Let- 
ters, 49(2):225-232, 1977. 

15 Daniel Sheppard, Penghao Xiao, William Chemelewski, 
Duane D. Johnson, and Graeme Henkelman. The Jour- 
nal of Chemical Physics, 136(7) :074103, 2012. 

16 J. P. Perdew, K. Burke, and M. Ernzerhof. Phys. Rev. 
Lett, 77:3865, 1996. 

17 Paolo Giannozzi, Stefano Baroni, Nicola Bonini, Matteo 
Calandra, Roberto Car, Carlo Cavazzoni, Davide Ceresoli, 



Guido L Chiarotti, Matteo Cococcioni, Ismaila Dabo, An- 
drea Dal Corso, Stefano de Gironcoli, Stefano Fabris, 
Guido Fratesi, Ralph Gebauer, Uwe Gerstmann, Chris- 
tos Gougoussis, Anton Kokalj, Michele Lazzeri, Layla 
Martin-Samos, Nicola Marzari, Francesco Mauri, Riccardo 
Mazzarello, Stefano Paolini, Alfredo Pasquarello, Lorenzo 
Paulatto, Carlo Sbraccia, Sandro Scandolo, Gabriele 
Sclauzero, Ari P Seitsonen, Alexander Smogunov, Paolo 
Umari, and Renata M Wentzcovitch. Journal of Physics: 
Condensed Matter, 21(39):395502 (19pp), 2009. 
Kyle J Caspersen and Emily A Carter. Proceedings of 
the National Academy of Sciences of the United States of 
America, 102(19) :6738-43, May 2005. 

B. S. Pujari, A. Tokarev, and D. Saraf. J. of. Phys. : 
Condensed Matter, 2012. In press. 

Sanghee Kim. Angewandte Chemie (International ed. in 
English), 48(42): 7740-7743, 2009. 

R. B. Heimann, J. Kleiman, and Salansky N. M. Nature, 
306:164-167, 1983. 

B. P. Stoicheff. Tetrahedron, 17:135-145, 1962. 

L. S. Nelson, A. G. Whittaker, and B. Tooper. High Temp. 

Sci., 4:445-477, 1972. 

B. V. Lebedev. Russian Chemical Bulletin, 49(6):965-975, 
2000. 

Stefan Grimme. Journal of Computational Chemistry, 
27(15):1787-1799, 2006. 

Vincenzo Barone, Maurizio Casarin, Daniel Forrer, 
Michele Pavone, Mauro Sambi, and Andrea Vittadini. 
Journal of Computational Chemistry, 30(6):934-939, 2009. 



